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Abstract: The effect of gravitational tidal forces on photon propagation in curved 
spacetime is investigated. It is found that the imaginary part of the local refractive 
index \m.n{u]uj) may be negative as well as positive, corresponding to a local am- 
plification as well as attenuation of the amplitude of the renormalized photon field. 
This is interpreted in terms of the effect of tidal forces on the virtual e^e~ cloud 
surrounding the bare photon field — a positive/negative Im?T,(M; u) corresponds to an 
increased dressing/undressing of the bare photon. Below threshold decays of the 
photon to e^e~ pairs can occur. Photon undressing in the vicinity of a black hole 
singularity is described as an example. These results are shown to be consistent with 
unitarity and the optical theorem in curved spacetime, which is derived here both in 
a local form and integrated over the photon trajectory. 



1. Introduction 

A consequence of vacuum polarization in QED in curved spacetime is that pho- 
tons do not simply propagate along classical null geodesic trajectories. Instead, the 
spacetime acts as a medium with a non-trivial refractive index. In a series of pa- 
pers [1-4], we have developed the theory of photon propagation in QED in curved 
spacetime and calculated the full frequency dependence of the refractive index in 
terms of a geometric quantity, the Van Vleck-Morette (VVM) determinant, which 
characterizes the null geodesic congruence around the classical trajectory. This has 
resolved the long-standing puzzle originating from the discovery of Drummond and 
Hathrell [5] that at low frequencies, the phase velocity may be supcrluminal. The 
geometry induces a novel analytic structure for the Green functions and refractive in- 
dex which modifies the conventional Kramers- Kronig dispersion relation and removes 
the apparent paradox between low- frequency superluminal motion and causality [3] . 

An initially surprising result of this work was that in some backgrounds, the 
imaginary part of the refractive index may be negative. In a conventional optical 
medium, Im n{u) is always positive, corresponding to scattering of photons from the 
beam and a reduction in the amplitude. A negative Im n{u) implies gain, with energy 
pumped in from an external source. In quantum field theory, Im n{(j) is related via 
the optical theorem to the rate of production of e'^e~ pairs and therefore vanishes; 
even for a massive "photon" above the pair production threshold, Im n{uj) would 
be manifestly positive. The fact that Im n{uj) can be negative in curved spacetime, 
and can be non-vanishing and positive below the e+e" threshold, requires a careful 
explanation and interpretation. In particular, it involves a reappraisal of the optical 
theorem in curved spacetime. 

In this paper, we summarize the resolution of this problem, leaving the full formal 
development to a longer paper [6] . We also provide an intuitive picture of how the 
amplitude may either increase or decrease depending on the sign of Im n{uj) based 
on the idea of "dressing" a renormalized quantum field. 

In QED, vacuum polarization means we may regard the renormalized photon 
field as the bare field dressed with a cloud of virtual e+e^ pairs. Since the electron 
has a mass, this cloud has a scale characterized by the electron Compton wavelength 
-^c = h/mc, so the renormalized photon should be thought of as an extended body 
with size Ac propagating through curved spacetime with scale L ~ where TZ 

is a typical curvature. Classically, massless point particles follow null geodesies in 
curved spacetime. Extended bodies, on the other hand, are subject to gravitational 
tidal forces because nearby geodesies can diverge or converge. In quantum field 
theory, because of vacuum polarization, we therefore expect that when L becomes 
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comparable to Ac, these tidal forces will distort the virtual cloud and affect the degree 
of dressing of the renormalized photon. 

The conceptual picture that we shall confirm is the following. To linear order 
in the curvature, when the rate of change of the acceleration of nearby geodesies 
is positive — in the direction of increased stretching — more virtual e^e~ pairs are 
produced and the photon's amplitude decreases. This is illustrated in the top diagram 
of Fig. m which indicates that the photon is becoming "more dressed" . If the rate of 
change of the acceleration is at least constant over a finite region along the geodesic 
then this can interpreted as real pair production and the amplitude reduces at a rate 
which is at least as great as exponential. On the contrary, when the rate of change 
of the acceleration of nearby geodesies is negative — in the direction of increased 
squeezing — virtual e~^e~ pairs are forced back into the photon. In this case the 
amplitude of is amplified. This is illustrated in the lower diagram of Fig. ^ 




Figure 1: Diagrams which illustrate the the cases where the photon amplitude decreases (top) 
and increases (bottom) arising from the production of more/less virtual e"'"e~ pairs, respectively. 
The top/bottom cases occurs when the rate of change of the acceleration on nearby geodesies is 
positive/negative. 

This mechanism gives rise to a perturbative contribution to Im n{u) of order 
ujdyTZ/m^ (where m is a lightcone coordinate along the null geodesic) in cases where 
the curvature is not constant along the photon trajectory. Im n{u) can be positive or 
negative and the photon amplitude can decrease (even below threshold) or increase. 
In many respects, the situation resembles the transient phase immediately following 
an interaction being turned on in conventional flat spacetime QFT. In both cases, the 
possibility of amplification as well as attenuation originates in the loss of translation 
invariance, which enables energy conservation to be circumvented. 

In addition to this very intuitive mechanism there is a further effect, non- 
perturbative in oj^/TZ/m?, which can give rise to an imaginary part for the refractive 
index. This is related to the occurrence of conjugate points along the photon's null 
geodesic 7. Two points on 7 which can be linked by a continuous family of geodesies 
arbitrarily close to 7 are said to be conjugate. The existence of these nearby geodesies 
implies there are zero modes in the fluctuations around the classical path and conse- 
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quently singularities in the VVM matrix [1-3] . This is the uhimate geometric origin 
of the novel analytic structure for the refractive index in curved spacetime and in 
some cases may imply the existence of a non- vanishing contribution to Im n{u}). This 
effect is present even in the simple example of symmetric plane waves, or Cahen- 
Wallach spaces [7] , where the curvature is constant along the photon path. In these 
cases, Im n{u) is always positive,^ a result which can now be understood from our 
analysis of the optical theorem in curved spacetime. For other backgrounds, there is 
no restriction on the sign of Im n{u). 

Clearly, the geometry of fluctuations around the classical photon trajectory is 
central to this analysis. Consider a null congruence centred on 7 and let z^{u), i = 1,2 
be the two spacelike connecting vectors linking neighbouring geodesies. The z'^{u) 
satisfy the Jacobi equation 

z' + hij{u)z^ = (1) 

where hij{u) = K^ujuiu), the components of the Riemann tensor in a frame associated 
to the Brinkmann coordinates (m, v, z*) adapted to the geodesic 7, which is the curve 
V = z^ = 0. The metric in the neighbourhood of 7 takes the form 

ds"^ = 2du dv - hij{u)z'z^ du^ + dz" dz^ . (2) 

This metric is known as the Penrose limit for the given background and null geodesic 
7 and encodes the geometry of geodesic deviation [8,9]. For later use, we also 
introduce the Rosen coordinates (m, V, x"), a = 1, 2, in which the approximate metric 
(H) takes the form 

ds"^ = 2du dV + Cab{u)dx'' dx^ . (3) 
The affine parameter u along 7 is common to both sets of coordinates. 



A key role is played by the Van Vleck-Morette determinant A(x, x') evaluated 
at the two points G 7 with affine parameters u and u' , so we write A{u,u'). 
There are different ways to define the VVM determinant but one very physical way 
is via solutions of the geodesic equation (p. If we write the solution for the Jacobi 
field z^{u) in terms of some initial data at u', 

z'{u) = ^'j{u,u')z^{u') + £/'j{u,u')z^{u') , (4) 

then [3,4] 

^ ^ det^(M,M') ^ ^ 



^Note that in refs. [3,4], there is an overall sign error in the refractive index for scalar QED. In 
particular, Imn(a;) for scalar QED in a Ricci-flat symmetric plane wave should be positive. 
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2. Vacuum Polarization and the Refractive Index 

The Penrose limit will be a good approximation for a general curved spacetime 
provided the curvature scale is greater than the electron Compton wavelength, i.e. 
L ^ Ac, or equivalently TZ/rn^ ^ 1, and we work in the geometric optics, or eikonal, 
approximation where the curvature scale is large compared to the photon wavelength, 
i.e. L » A, or TZ/u"^ <C 1. We refer to these two conditions as the "WKB limit". 
This leaves the dimensionless combination uj\/TZ/m^ of the three scales as a free 
parameter in our calculation of the full frequency dependence of n{uj). 

To simplify the presentation and focus on the key issues of the sign of Im n{uj) 
and the optical theorem, we quote results here for a model field theory compris- 
ing a massless scalar (photon) A interacting with two massive scalars (electrons) 
through an eAcf)^ coupling. The case of QED involves some further subtleties with 
renormalization and is described elsewhere [6]. 

In the leading eikonal approximation, the solution of the free wave equation 
□ [x] — takes the approximate form 



where g = det^^^j, and k'^ = d^'d{x) is the tangent vector to a null congruence 
g{k, k) = 0. Along the geodesies '(9(a;) = V" is constant while the affine parameter u 
varies. If we pick out 7 to be the geodesic with V = then in the neighbourhood of 
7, V is identified with the Rosen coordinate. 

In our model field theory, vacuum polarization is given by the usual QED-like 
one-loop graph of C(e^) with the scalars in the loop and external A lines. Since 
the usual notion of a particle is inherently ambiguous in a general curved spacetime, 
we shall focus on the behaviour of the field itself and use the techniques of non- 
equilibrium QFT to set-up an initial value problem. We therefore imagine that on 
the null surface u = Uq the coupling constant e is turned on instantaneously and 
then watch as the bare field evolves along the null direction u.^ The field satisfies a 
non-local equation of the form 



where IIr{x,x') is the one-loop retarded self-energy. Since the latter is only non- 
vanishing when x' lies in the past light-cone of x, the equation is manifestly causal. 

^Strictly speaking a null surface is not a Cauchy surface. However, in the philosophy of light- 
cone quantization one can imagine defining boundary conditions in the transverse directions {V, x"") 
which have the eflfect of giving a well-defined initial value problem. 



(6) 
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The aim is to compute the solution of this perturbatively in e starting with the 
classical solution (P) in the form 

= A(°)(x)(l + ze2Q(i)(x) + ---) . (8) 

The one-loop retarded self-energy is 

Ur{x, x') = i9{u - u') Im Gf{x, x'f , (9) 

where Gf{x.,x') is the Feynman propagator of and in order to extract the 

one-loop correction we need to calculate the integral 



(iV V^Mni?(x,x')A(°)(a;') • (10) 

This integral becomes tractable in the WKB limit, where we can replace a general 
background spacetime by its Penrose limit (H) or (j^). 

First, we introduce some notation which will be useful when we discuss the 
optical theorem. The general on-shell modes for a classical field of mass m in a plane 
wave background in Rosen coordinates are 

0p(x) = (7(M)-i/^exp [luV + ipax'' - ^{r\u)paPb + m\)] , (11) 

where p = {u),pa) includes a transverse momentum and we have defined the 2x2 
matrix 

iIj{u) = / du'C{u')-^ . (12) 



This matrix is related to the VVM determinant via 

1 (u - u'f 

""'^ ^ ' detain, u') ■ (^2) 

These modes satisfy the orthonormality property, 

(0p,, (j)p)=tj d?x^ 4>;,{x)K<Pp{x) = 2uj{2Txf6{uj - uj')6^^\pa - p'a) , (14) 

where the integral is over x = {V,x'^). The original massless photon modes A^^^ of 
(^ were chosen with Pa = 0. 

Now define the partial "Fourier transform" of the self-energy with respect to 
these (massless) modes: 

{2-k)'6{u - {pa - p'a) flRiu, U'- U, Pa) 

= j £xy/g j £x'^' (l)l{x)TlR{x,x')(l)p,{x') . ^^^^ 
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Solving the field equation and taking into account that IIr{x,x') has support 
only for u' < u, we can then write the quantum correction Q^^'^ of (Q) in the form: 

Q^'\u) = — / du" / du'IlR{u'\u';u,0) . (16) 



Returning to ([T0|), we now summarize some of the key points of the calculation, 



while the full details can be found in [6]: 

(a) In the WKB limit, for a fixed x on the null geodesic 7, the integral over x' is 
dominated by a saddle point with x' G 7 in the causal past of x. To leading order, 
the metric can be replaced by the approximate plane wave metric @) or (|]) which 
captures the fiuctuations around 7 to Gaussian order. 

(b) The propagator for a massive scalar field in a plane wave is known exactly [10]: 

G^{x,x')=^y^A(^)j^ J^2^~"^^^- (17) 

where a{x,x') is the geodesic interval between x and x' and the saddle point over 
the x', T\ and T2 (one for each propagator) integrals is at 

T T 



V' = x"' = , u' = u- 2m^u 



T1+T2 



What remains is an integral over the Feynman parameter < ^ = j^jt < 1 and 
the point x' on 7 labelled by affine parameter u' with Uq < u' < u. Performing 
the Gaussian integrals around the saddle-point gives the equation for the quantum 
correction Q^^\x). Evaluating this along 7 gives 



1 r \/Mu U') rm^(u-u') 

which can easily be integrated with initial condition Q^^\uq) = giving 

Q^'\u) = -V / du" / du' / d^e-^^ . (20) 

The curvature is encoded in the Van Vleck-Morette determinant. The result (|19]) is 
divergent for two entirely different reasons: (i) there is the standard UV divergence 
as u' — > u. This is the same divergence as in fiat space that is absorbed by a mass 
counterterm for the A field. This divergence does not affect the the imaginary part 
of Q^^^ which is our main focus here, so we will not explain this point any further; 
(ii) there is a divergence whenever u' and u are conjugate points coming from the 
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VVM determinant. These divergences are treated by implementing the real space 
Feynman prescription, u — u' ^ u — u' — iO'^. 



Also note that because of the symmetry A{x,x') = A{x',x), the imaginary part 
of Q^^^ (u) is given simply by extending the range of the integration from Uq < u' < u" 
to Uq <u' < u. We find 



ImQ(^)| 



u] 



du" / du 



A(u",u' 



u — u 



die 



(21) 



'UQ JUQ "'Jo 

This quantity appears in the optical theorem, as described in section 3. 



(c) In general the one-loop correction can grow as a function of u and it becomes 
necessary to perform a re-summation using the dynamical renormalization group [11]. 
This results in an exponentiation of the result and the solution to the quantum 
corrected field equation (^ is then 



Aiu) = ^7(^z)-i/V'^^e^'='s(^)H 



(22) 



The result can be interpreted in terms of a position-dependent refractive index 

n{u;uj) = l + -dM^'\u) (23) 

CO 

with a real part that describes a non-trivial correction to the phase velocity of the 
massless field A{x) and an imaginary part that describes a non-trivial correction to 
the amplitude 



\A{u) \ = \A{uo)\exp 



u 

—Lu I du lmn{u; u) 

UQ 



(24) 



along 7. 



The result (^) is already non-trivial even in flat space. This describes a tran- 
sient region where the bare field A a.t u = uq renormalizes and becomes a dressed 
field in real time. Indeed, as u increases, lmduQ^^\u) settles down after some oscilla- 
tions to a positive constant and the ratio |y4(u)/y4(Mo)| approaches the wavefunction 
renormalization factor that one would calculate in the equilibrium theory: 



Aiu) 



A{uo) 



—J- exp 



1 - 



+ ■■ 



(25) 



Notice that in A(j)'^ theory in d = 4, wavefunction renormalization is a finite effect. 
The decrease in the renormalized amplitude can be interpreted in terms of the in- 
tuitive picture of the bare field surrounded by a virtual cloud of pairs. As this 
dressing increases, the amplitude falls. 
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We can decouple this transient initial dressing of the field from the effective of 
gravitational tidal forces by taking uq — ?■ —oo and then turning on the curvature at 
some finite u after the transients have died out. However, we must not forget the 
initial dressing process which is essential in maintaining consistency with unitarity, 
as explained later in our discussion of the optical theorem. 

In the weak curvature limit 7^ ^ m, we can expand the VVM determinant to 
linear order in the curvature 

A(u, u') = l + lRuu{u)iu - u'f - ^duRuuiu)iu - u'f + ■■■ (26) 

where the Ricci tensor Ruu = R\iu- Substituting this expansion into the integral 
(|T9|), and after UV regularization, we determine the correction to the refractive index 

n{u; u) = l- ^ Ruu{u) - ^ duRuuiu) , (27) 

In this approximation, the quantum correction to the amplitude relative to a point 
u' ^ Uq is 



\A{u)\ = \Aiu')\exp 



4207r2m6 



{Ruu{u) - R 



(28) 



If the Ricci curvature component Ruu{u) decreases along 7, then from the Jacobi 
equation (p we see that the rate of acceleration of the nearby geodesies is positive, 
i.e the tidal forces are increasing in the direction of stretching the virtual cloud. In 
this case, the amplitude decreases because more virtual pairs are being produced: 
A is becoming more dressed. On the contrary, if Ruu{u) increases along 7, i.e the 
tidal forces are increasing in the direction of squeezing, then the amplitude increases 
because virtual pairs are interacting to produce more A: A is being "undressed" and 
is returning to its bare state. These effects are illustrated heuristically in Fig. |l[ 



We now come to the second, non-perturbative, mechanism which can produce 
an imaginary part for the refractive index. If Ruu{u) is constant, in which case 
the approximate metric (|]) and (^) is a symmetric plane wave, or Cahen-Wallach 
space [7], then to linear order in the curvature the amplitude is constant. In fact, 
Imn(a;) = to all orders in the curvature expansion. However, in the case where 
at least one of the eigenvalues of the constant hij is negative there is a positive 
contribution to lmn{u), so a decaying amplitude, which is non-perturbative in the 



curvature. The simplest case to consider is hij 
determinant is 



a^6ij. In this case, the VVM 



A(m, u) 



sinh a(u 



u 



u' 



(29) 
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Introducing t = u — u' we have 



n{u) = l + -——J ——- d^e-^^ . (30) 



327r^Ci;^ 7o sinh(crt) 

lmn{u) can be computed by rotating the t = u — u' integral t — —it + iO"'' and 
picking up contributions from the poles att = — ,n = l,2,.... This gives 



-,2 /■! 



Imn(a;) = — ^/ d^ ,,,2 , (31) 

which is, indeed, non-perturbative in the curvature. In this case the result is inde- 
pendent of u and so it implies a constant rate of production of pairs. Such a process 
is kinematically disallowed in fiat space, but there is no such threshold constraint in 
curved space. 

This result leads to a very non-trivial check of the formalism in the following 
way [6]. If we consider the space dS^ x R with metric 

ds"^ = -dt^ + cosh2(at) (rf^^ + sin^ 9 dcj)'^) + dx^ , (32) 

then the plane wave limit about the geodesic t = x and a fixed point on the S*^, 
e.g. ^ = ^, = 0, is obtained by transforming 

and then expanding in powers of v and = 6* — |, = 0: 

ds^ 2du dv + cosh^ {^)dx^ dx" . (34) 

This is a symmetric plane wave in Rosen coordinates (|^) with hij = The 
non-trivial check on (^) comes about because the decays of scalar fields A (jxp 
can be calculated exactly in de Sitter space due to large amount of symmetry as a 
tree-level process for any masses. In our case A is massless and using the results 
in [12, 13] and taking the limits u ^ a and m ^ a one gets a decay rate which is 
precisely 2c<;Imn(a;) as found in (^Tj). 



We therefore see that in curved spacetime, the effect of gravitational tidal forces 
on the virtual cloud surrounding the bare field means that the photon can become 
increasingly "undressed" as well as dressed as it propagates through spacetime. An 
especially interesting example of this gravitational undressing phenomenon is photon 
propagation near a singularity, for example a black hole. The Penrose limit near the 
singularity has the universal form of a homogeneous plane wave [9] 

h.A^) = ^-^k, , (35) 
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with Ui being either real or imaginary (and choosing a, > in the former). For 
example, the near-singularity Penrose limit of a null geodesic with non-vanishing 
angular momentum in the Schwarzschild solution is described by such a plane wave 
with «! = I and 02 = |. As the singularity is approached m — )■ (with m < 0) the 
u' integral is dominated by its endpoint u' = uq = —00 and we find 

ImQW(M) — ^ -\u\-^^ . (36) 

The amplitude therefore increases and the photon becomes progressively undressed as 
the singularity is approached. Notice that the amplitude remains finite even though 
\m.n{u\ uj) itself diverges if ai + 02 < 2 as in the case of Schwarzschild. When either 
of the Ui is imaginary the amplitude has a oscillating behaviour. 



3. Optical Theorem 



The final issue we need to address is how the optical theorem and unitarity can 
be satisfied if the amplitude can increase and Im n{u, u) can be negative. The answer 
to this puzzle is that we have only shown that the amplitude can increase locally in 
a region well away from the initial-value surface. In order to verify that unitarity 
is satisfied, and that the decrease in the amplitude is related to the probability of 
the production of pairs, we have to consider the amplitude over the whole region 
from the initial-value surface uq to u. It turns out that while the amplitude can 
increase locally, the overall amplitude as measured from uq always decreases. The 
initial transient dressing of the A field plays an important role since it is only a field 
that is already dressed that can increase its amplitude by becoming more undressed; 
a bare field, however, cannot be undressed any further! 

To derive the optical theorem, note first that 

ru g2 ru pu" ^ ^2 

/ du" lmn{u";uj) = — - / du" / du lmllii{u" ,u'; u,0) = — lmQ^^\u) . 
Juo Juo Juo ^ 

(37) 

where Im Q^^\u) is given in eg . (|2l|) . Now consider the the transition probability for 
the tree-level process A — (pcj): 

d^Pi d^p2 



■ (38) 

Here, 0p(a;) are the massive on-shell modes defined in (pTj) and the integrals are 
over the three-dimensional p = (uj,pa)- The normalization factor J\f = {A^'^\ A^^^) = 
2u6^^\0) cancels against an overall integration on {V, x"). Recall that the Wightman 
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functions G+{x,x') = {0\(f){x)(f){x')\0) and G-{x,x') = {0\(j){x')(j){x)\0) are given in 
terms of these modes by 

G^{x,x') = ± [ <p^{x)<p,{xy . (39) 

and that the frequencies are all positive since they refer to physical particles in (|38|) . 
The retarded self-energy (^ is expressed in terms of the Wightman functions by 

Ur{x,x') = ie{u-u')[G+{x,x'f -G4x,x'y] . (40) 

A short calculation using the properties G+(x,x')* = G-{x,x') and G'+(x',x) = 
G-{x,x'), and noting that J dV'G^{x, x'^A^'^^x') = for positive u, now shows 
that 

^2 /•« 



Pa^4,4,{u) = — du" / du' lmUR{u" ,u;u,0) , (41) 

^ Juo Juo 



and so 



PA-^<t,<t}{u) = 2u du" lmn{u";u) . (42) 

Juo 

This is the statement of the optical theorem in curved spacetime. Unitarity is re- 
spected and the integral over the whole trajectory of lmn{u;u) is positive. This 
shows how unitarity prevents a local amplification of the amplitude from becoming 
unbounded. 

We can also write a local version, defining T(u) = duPA~¥4>4,{u) as the "instanta- 
neous rate" of pair production. Then, 

T{u) = 2u}lmn{u;uj) = — / du lmIlR{u,u';uj,0) . (43) 



Unlike the integrated form, there is no positivity constraint on ( ^3D in general. This 
explains why there is no confiict with unitarity in the examples where we have found 
lmn{u]u) < 0. However, in cases where we have translation invariance along the 
photon trajectory, as in fiat spacetime, n{u) is u-independent — as long as one is 
beyond the transient region — and Pa-s>(/>(/)(m) is proportional to u, so unitarity implies 
Im?T,(a;) > and F > can be properly interpreted as the rate of A — )■ (fxp. 



4. Conclusions 



In conclusion, we have shown how the effect of gravitational tidal forces on 
vacuum polarization can alter the dressing of a photon as it propagates through 



space. In particular, we have seen how the imaginary part Im?T,(n; u) of the position- 
dependent refractive index can be negative as well as positive, corresponding to 
"undressing" of the photon rather than the conventional dressing. Two mechanisms 
were identified. The first, of order uodyTZ/m^, admits an intuitive interpretation 
in terms of curvature variations along the photon's trajectory altering the balance 
of the bare field with its virtual e^e~ cloud, with increasing stretching (squeezing) 
giving rise to more dressing (undressing). The second, which is non-perturbative in 
uTZ/m?', can occur even when the curvature is constant and is related to the existence 
of conjugate points on the photon's null geodesic. 

Nevertheless, unitarity is still respected and the optical theorem still holds in 
curved spacetime. In its integrated form (|4l|) , (|4^) , it relates the total probability 
for e"'"e~ pair production to the integral of lYsin{u;uj) along the photon trajectory, 
which is manifestly positive. Except in the special case of translation invariance along 
the null geodesic, the corresponding local form (^) has no positivity constraint and 
the usual interpretation of Imn(M; u) as the rate of pair production can break down. 
However, it does describe the variation of the amplitude and its interpretation in 
terms of dressing and undressing of the photon field. In a sense, photon propagation 
in curved spacetime resembles the initial transient phase in fiat spacetime, with the 
characteristic features of an oscillating amplitude and below-threshold decay. 

In this paper, we have described our results in terms of a scalar Acj)^ model of 
QED. In QED itself, there are two further complications — the photon has a polar- 
ization, whose direction influences the dynamical dressing, and there are subtleties 
related to renormalization, requiring a reformulation of the initial-value problem to 
circumvent the short- distance divergences at the initial value surface caused by the 
lack of asymptotic freedom of QED. These issues are explored in a more compre- 
hensive treatment of the effect of gravitational tidal forces on renormalized quantum 
fields in ref. [6]. 
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